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Abstract In this work we extend previous work on the evo- 
lution of a Primordial Black Hole (PBH) to address the pres- 
ence of a dark energy component with a super-negative equa- 
tion of state as a background, investigating the competition 
between the radiation accretion, the Hawking evaporation 
and the phantom accretion, the latter two causing a decrease 
on black hole mass. It is found that there is an instant during 
the matter-dominated era after which the radiation accretion 
becomes negligible compared to the phantom accretion. The 
Hawking evaporation may become important again depend- 
ing on a mass threshold. The evaporation of PBHs is quite 
modified at late times by these effects, but only if the Gen- 
eralized Second Law of thermodynamics is violated. 

Keywords Black holes; Cosmology; Hawking radiation; 
phantom energy. 

PACS 04.70.-s; 98.80.-k; 95.35.-^d 



1 Introduction 

The now widely accepted accelerated expansion of the Uni- 
verse in its recent history is yet to be fully explained. Several 
possibilities to reproduce this effect have been advanced, 
ranging from "conservative" to very unusual ones requiring 
new physics. One of the most economical hypotheses that 
has received a great deal of attention is the late dominance 
of a fluid with an "anomalous" equation of state, a sort of 
analogue of the inflationary proposals but at a lower energy 
scale, the so-called dark energy. As is customary to write a 
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generical equation of state in the form P = cop and in spite 
that values of (0 larger than — 1 are usually considered, some 
works have raised the possibility that the dark sector may 
be characterized by a fluid with an equation of state with 
ft) < — 1 , known throughout the literature as the phantom 
energy. 

There are many physical consequences of such phantom 
component in a variety of physical species present in the 
Universe, most notably the spacelike singularity known as 
the Big Rip LL^, or even more fabulous possibilities, like 
the Big Trip |I3]|4]|. Some effort has been made to remove the 
Big Rip singularity, but it is still premature to rule out or 
support definitely any scenario. 

We work within a general phantom energy scenario in 
this paper. It has already been acknowledged that, being such 
an exotic physical species, the phantom energy may also 
change the accretion regime of black holes 1 5 1 . In the present 
paper we investigate the influence of phantom energy accre- 
tion onto primordial black holes (hereafter PBHs) together 
with the radiation and matter accretion/evaporation formerly 
addressed. 

The PBH interaction with different types of energy in 
the universe is the continuous subject of several sudies, as 
well as their interaction with cosmological boundary condid- 
iotns [6|. Several numerical results also work as test fields 
for altemate gravitational theories, and the questions regard- 
ing their very formation at extreme cosmological scenarios 
are beginning to yield several interesting results |7|. 

We shall focus on the new features specifically intro- 
duced by the phantom era [8|, and generally refer to the 
full evolution of the PBHs across the mass-time plane. Pre- 
vious attempts to address this problem have been limited 
to the consideration of the black holes plus phantom fluid 
only, although there is a more subtle interplay among com- 
ponents when radiation and matter are also included (as will 
be shown below). It is also of interest to revisit the issue 
of the black hole behavior in the radiation-dominated and 
matter-dominated eras (that is, well before the phantom com- 
ponent can be important) for a complete assessment of the 
fate of PBHs, especially their behavior in the matter-dominated 
era where dark matter may fuel their growth. 



2 PBHs evolution in tlie early Universe 



2.3 Accretion regimes 



2.1 The radiation-only accretion equation 

Our starting point will be the evolution equation for PBHs 
in the radiation-dominated era addressed by several authors 
(see Custodio and Horvath |9 1 and references therein), which 
takes into account the accretion of radiation and the Hawk- 
ing evaporation at a semiclassical level. Ignoring the (poten- 
tially relevant) "grey factors" in the absorption of radiation, 
the resulting differential equation for the black hole mass M 
reads quite generally 



dM 
dt 



A{M) 21nG^ 



Pr.6iT)M^ 



(1) 



with t being the cosmological time, A{M) = ^a{M), with 
a{M) called the running constant IIIOII . counting the degrees 
of freedom of the emitted particles on the Hawking radiation 
(in CGS units, A = 7, 8 x 10^^ g^/s for black holes evaporat- 
ing today nil), and PndiT) the radiation energy density at 
temperature T at the time f . 

In a Universe also filled with phantom energy, the accre- 
tion of such exotic component should also be taken into ac- 
count. Babichev, Dokuchaev and Eroshenko |5 1 have worked 
out a differential equation for a black hole accreting phan- 
tom energy only, obtaining a counterintuitive result that phan- 
tom energy accretion decreases the overall black hole mass. 
The expression is similar to the accretion term in eq. ([T|, and 
is given by 



dM 
dt 
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-P{pph) 



(2) 



2.2 The complete accretion equation 



Considering the radiation accretion and evaporation terms 
from eq. ([TJ together with the new phantom energy accretion 
term in eq. (O, and assuming no interaction between the two 
different species, the complete equation for the accretion of 
the different types of energy into the black hole is just 



dM 
dt 



-^^ + ^ [277rft,d(r) + 167r (pph +/?(Pph))] M^ 

(3) 

Using for the phantom energy p{p) = cop, (O < —1, the 
phantom component of the accretion may be written as 



Pph+P{pph) = (l+«)Pph 
and the complete accretion equation becomes 



(4) 



dM 
dt 



A{M) G^ 



[277rp,ad(r) + 1671(1 + ft))Pph] m2 



As is well-known, the Friedmann equation can be solved to 
follow the cosmological evolution of the phantom energy, as 
given by Babichev, Dokuchaev and Eroshenko 13 . 



\p+p\ °^a 
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(6) 



neglecting all other contributions. The densities of the 
radiation, matter and phantom energy terms evolve as repre- 
sented in the graphic shown in Figure [TJ 



log(p,) 




log(a) 



Fig. 1 Evolution of the radiation, matter and phantom energy densities 
with the scale factor. The exact position of fpj, depends on the densities 
and equation of state of the radiation and phantom energy. 



As expected, there is an epoch in which the radiation and 
phantom energy accretion terms from eq. (|5]l become com- 
parable. We call such an epoch the phantom time, or fph. It 
must be noted that this instant is distinct from the one when 
the lines of Figure [T] cross each other. The phantom time 
represents the cosmological instant when the phantom en- 
ergy accretion term dominates the radiation term, changing 
drastically the black hole evolution dynamics. We can calcu- 
late the value of this time as a function of the initial radiation 
and phantom energy densities. 

The radiation density as a function of the scale factor is 
given by the Friedmann equation, piaj = p^^^ (22.^ During 
the matter-dominated era, the scale factor as a function of 
time is given by 



a{t) 
ao 



m"' 



(7) 



Therefore, the radiation density evolves in the matter- 
dominated era as 



Prad — Prad 



m 



(8) 



(5) 



Similarly, with the phantom energy eq. ^ and evolving 
according to eq. (|6]l, and with the time dependence of the 



scale factor evolving as of eq. (|7]l, the phantom energy den- 
sity as a function of time is 



Pph 



Pp°h 



1- 
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(9) 



The epoch when the phantom energy accretion is as im- 
portant as the radiation accretion is the instant when, equat- 
ing both expressions according to eq. (|5]l 



16,, 

Prad = -^(l + «)Pph 



(10) 



Inserting the time dependences calculated in eq. (O and 
eq. (|9]l, this equation yields ths phantom time. 



(1+Z)dec 



m 



(14) 



Well before the phantom energy becomes important, the 
PBH mass equation, including now the dark matter contri- 
bution, is just 



dM 
IT 



A 



^— PradM - 
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= — - — M 



(15) 



We must remark that we are always referring to a diffuse 
CDM component, an appropriate assumption prior to any 
structure formation. 
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with Hq expressed in ^-^^ and, as the initial values p^^ and 

Prad' calculated at the end of the matter-dominated era. 

We can express this transition time in terms of the red- 
shift, using eq. jlQI) . with the initial conditions p^^^ = 8. 12 x 



10-13 £Ig and p° = 1.79 x 10"^ -^ appropriate for the 
obtained conditions, finally coming to Zph ~ 3.1. 

It is reasonable to suppose the transition between radia- 
tion and phantom accretion to be instantaneous due to the 
very steep radiation/phantom density ratio, which can be 
easily seen by rewriting eq. jilt for an arbitrary epoch. 



3.2 Numerical predictions 

Because we are interested in the fate of a wide range of black 
hole masses, we should integrate equation ( 115b numerically 
for several initial conditions and cosmological parameters. 

To solve this equation, we first rewrite it in explicitly 
time-dependent terms 



Pm = Pm(l+z) 



Pde 



(l+^)dec 



;i+z)'=Pdec(^)' (16) 



Letting m = 100 GeV and Tdec ^ 1 GeV, (1 +z)dec 
4.26 X 10'2 yields 



Pl^ = E^U+co\('J^) 



h2(l+m) 



Pph Pph 



V 2 / 



(12) 



The radiation density quickly becomes negligible com- 
pared to the phantom energy. The higher the \co\, the quicker 
the transition becomes. 
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Setting pO = ^^n, = 
■, equation ( 115b reads 



32nGt 



2.24 X 10 



-30 



(17) 



and Piad = 



3 Effects of dark matter accretion 

3.1 General results 

Up to this point we have neglected completely the possible 
effects of (cold) dark matter on the PBHs, which is a popular 
and reasonable explanation for the structure formation prob- 
lem. Within the CDM scenario, right after the decoupling of 
dark matter its accretion onto black holes will depend on the 
black hole cross-section for point-like particles. Therefore, 
the time dependence of the mass would be given by 



dM 
dt 



\6nG^ pm, n 



(13) 



where u^ is the dark matter particle density, computed after 
the decoupling 



dM 
dt 



A 
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(18) 

M 



32 c3 f2 c3 0.173c 

which can be solved introducing new variables y = ■§-, 

Mq = ^^ and X = log ( /- ) , yielding an equation of the 
form 



"^J -2 X 



-a2y^e' 



aj,y 



(19) 



with a\ = 



AG 



1.30xlO~ 



, 02 = i\a = 2.53125a and 



fl3 = 1.38x10-42^0. 

The dark matter accretion should be taken into account 
for X > jCdec- We may also introduce an instant x* similar to 
the phantom time, in which the dark matter and radiation 
accretion have the same value. An estimate for jcjec and x* is 
given by 



fo = 2.47 X IQ-^^Mo -^ Xdec = log ( ^'^^."^ ^°^ ) (20) 



ajy^e ^ = aj,y^ — > x* = log 



Mo 
1.83x1041 
Mo 



(21) 



Table [T| summarizes a few numerical estimates for x^g_, 
and x^,, covering most of the important PBH masses. 



Table 1 Numerical values for xjec and xq for some initial values for the 
black hole mass, along with calculations for the times of evaporation 
with (x^y^p) and without (xevap) dark matter 



Mo(g) 


-tdec 


X, 


-^evap 


m 

-^evap 


10« 


52.56 


76.59 


63.99 


63.99 


10" 


50.26 


74.29 


68.59 


68.59 


lO'o 


47.96 


71.98 


73.20 


73.10 


10" 


45.65 


69.68 


77.81 


77.81 


10'2 


43.35 


67.38 


82.41 


82.41 


10'3 


41.05 


65.08 


87.01 


87.01 


10" 


38.75 


62.77 


91.62 


91.62 


1015 


36.44 


60.47 


96.23 


96.22 


1016 


34.14 


58.17 


100.83 


100.83 



An inspection of Table \T\ shows that only black holes 
with masses greater than 10^ g should be influenced by the 
dark matter accretion at early times. However, this effect 
of the dark matter term happens to be small, because it is 
rapidly overcome by the accretion of radiation. This can be 
expected on physical grounds because the geometrical dilu- 
tion of the dark matter component "starves" the PBHs by 
quickly diminishing the flux of particles coming into them. 
Note that this particular evolution does not refer to much 
later epochs where dark matter halos had formed, possibly 
then contributing to the growth of PBHs as seeds for the ul- 
timate supermassive galactic residents. 

The numerical results for the evolution through time are 
depicted in Figure |2] for the highest initial condition, as an 
example. The resulting bump in the mass (Fig. |2]l has been 
exaggerated for the sake of clarity. 



4 Behavior of the critical mass function 

With expression eq. (lilt for the time, we can calculate the 
value of the critical mass Mc in the instant fph- From Custodio 
and Horvath 19], the expression for the critical mass is 



M,(0 - lOMHa 



1) 

Is/ 



(22) 



During the late phantom energy accretion dominance era, 
a critical mass function would be meaningless, since there 
is no longer a relevant mass increase mechanism. Thus, the 
largest value reachable by the critical mass in a Universe 
filled only by radiation and phantom energy is 



M" 



lOMHa 



3//n 



Pad , 



2(1+(U) 



(23) 



After this time, the Hawking evaporation is no longer a 
relevant mechanism for black hole mass decrease, until its 
mass reaches the transition value discussed in section|5] 

It is also convenient to calculate the initial mass of the 
black hole which disappears at fph. For that purpose, it is 
enough to consider only the Hawking term in eq. ([Tl, which 
yields the well-known solution 



1 



3A(M) 



Mi 



(24) 



where T is the evaporation timescale. Restoring the cgs units 
T reads 



10 



71 



Mi 

Mq 



(25) 



Combining eq. ( l22l) and eq. (|24] |. we find a third degree 
equation in Mc, whose solution is the critical mass of the 
black hole that will evaporate completely at f = fph 



Mi 



Mi 



3A{M) 



ioom2^^ 



fph 



(26) 



We use the numerical values of A(M) < 7,8 x lO^*" \ S 
and Mhsw = lOi^ g, as well as the numerical values of Pph, 
Prad, (0 and Hq necessary to compute fph- The instant when 
the critical mass assumes this value is found by inverting 
eq. (|22]|. 

Since the mass gain due to radiation accretion is not sub- 
stantial |9|, all black holes with M, < M^ , which reach crit- 
ical mass at Tcross ^ tc, will disappear before fph and will 
never reach the phantom era. 



5 The competition between phantom accretion and 
Hawking evaporation 

We have emphasized before that, since after fph there is no 
efficient mechanism that could increase the mass of the black 
holes, there is no longer a critical mass function. However, 
due to the presence of a phantom field, there are now two dis- 
tinct regimes of mass decrease, whose relative importance 
depends on the mass of a given PBH entering the phantom 
era. 

Taking eq. ^ and neglecting the radiation term, we can 
describe the evolution of black hole masses during the phan- 
tom era. Let us define a ratio between the two remaining 
terms. 



^(M): 



Mph _G M67r(l + ft))pph 
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M^ 



(27) 



or, in terms of a transition mass 



^{M) = 



with 



Mt = 



M 



A{M) 



l6nG^ (l + a))Pph 



1/4 



(28) 



(29) 



Substituting numerical values for the constants, we ob- 
tain an expression for M, in terms of the phantom field den- 
sity 



M,^5.5xlO"[(l + ft))Pph]-'/*g 



(30) 



with Pph given in g/cm-'. 

Since both regimes are of mass decrease, the black hole 
mass will diminish mostly due to phantom accretion un- 
til it reaches M,. After this, the predominant effect will be 
Hawking evaporation, since eq. ( 1281) shows that the change 
in regimes is sufficiently sudden for us to make this approx- 
imation. 

To find the time dependence of the transition mass, we 
must first know the evolution of the phantom density. Ac- 
cording to the Friedmann equations for the phantom fluid, 
we finally obtain Q 



(pph) 
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+ co) /8kG\ 
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(31) 



with (Ppjj)~^ being the initial density of the phantom field. 
Inserting this result on equation eq. ( |30l l the time depen- 
dence of Mt is obtained 



M,^ 



8.29 X 10^1 

(1 + 0))^ 



(pS 



ph^ 



3(1 + ©) /87rG 



3 



(32) 

The initial value of the transition mass depends on both 
the initial value of the phantom density and on CO. It is worth 
remarking that this transition mass is meaningless in the 
radiation-accretion regime. 

The differences between the three regimes is depicted in 
Figure 121 

It is important to stress that the Hawking evaporation 
does not become negligible after fph if taken into account 
as an independent process. However, the masses for which 
it becomes important (M <M,) drop by a factor of 10^ after 
the transition. This suddenly drives many black holes, but 
not all, into the new regime. When, however, the black holes 
reach the Planck mass, a full quantum gravity analysis be- 
comes necessary to properly determine its fate, since it has 
been shown that the Hawking evaporation no longer behaves 
as expected on such scales II13II14II . 




Fig. 2 Primordial black hole evolution in the matter-radiation- 
phantom energy scenario. The thick lines represent the different trajec- 
tories of black holes of different initial masses. The Big Rip singularity 
occurs at tg^. 



6 Conclusions 

We have studied in the present work the evolution of PBH 
for various regimes of accretion/evaporation in the very early 
and contemporary Universe. In particular, we have extended 
and clarified the evolution in the radiation-dominated and 
matter-dominated eras, including the features of diffuse CDM 
accretion producing only a small bump in the mass of the 
PBHs at early times. We have generally confirmed previ- 
ously known features of the semiclassical pictures of PBH 
evolution from a general point of view. Novel features are 
introduced in this scenario when a phantom energy compo- 
nent is introduced, as suggested by Babichev, Dokuchaev 
and Eroshenko 151 . 

Broadly speaking, a phantom field introduces another 
evaporation regime that competes with the celebrated Hawk- 
ing evaporation. We have found that the joint consideration 
of the relevant terms quenches the asymptotic approach to a 
common mass resulting from the phantom term only. This 
conclusion should, however, not be considered as definitive. 
Its validity rests on the assumption of the entropy for the 
phantom fluid being negligible, which is not the most gen- 
eral possibility. In fact, the enforcement of the Generalized 
Second Law (GSL) of thermodynamics would forbid the 
evaporation of the PBHs by phantom accretion |15 18] In 
addition, it is not clear whether the GSL should be valid in 
presence of the phantom fluid not respecting the dominant 
energy condition, as pointed out by Izquierdo and Pavon 
[15] , and models may be constructed in which the GSL must 
be modified. There is a rich variety of behaviors II16II17II 
within phantom energy models that remains to be explored 
in connection with the PBH evolution problem. In particular, 
late evaporation may conflict with the generalized second 
law of thermodynamics lil8J . 
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